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ABSTRACT 

Kinematic surveys of the dwarf spheroidal (dSph) satelHtes of the Milky Way are revealing 
tantalising hints about the structure of dark matter (DM) haloes at the low-mass end of the 
galaxy luminosity function. At the bright end, modelling of spiral galaxies has shown that 
their rotation curves are consistent with the hypothesis of a Universal Rotation Curve whose 
shape is supported by a cored dark matter halo. In this paper, we investigate whether the 
internal kinematics of the Milky Way dSphs are consistent with the particular cored DM 
distributions which reproduce the properties of spiral galaxies. Although the DM densities 
in dSphs are typically almost two orders of magnitude higher than those found in (larger) 
disk systems, we find consistency between dSph kinematics and Burkert DM haloes whose 
core radii ro and central densities po lie on the extrapolation of the scaling law seen in spiral 
galaxies: log po - « log ro + const with 0.9 < a < 1.1. We similarly find that the dSph data 
are consistent with the relation between po and baryon scale length seen in spiral galaxies. 
While the origin of these scaling relations is unclear, the finding that a single DM halo profile 
is consistent with kinematic data in galaxies of widely varying size, luminosity and Hubble 
Type is important for our understanding of observed galaxies and must be accounted for in 
models of galaxy formation. 

Key words: dark matter — galaxies; dwarf spheroidal — galaxies; kinematics and dynamics — 
Local Group — stellar dynamics 



1 INTRODUCTION 

In the current cosmological paradigm, dark matter (DM) 
provides the gravitational potential wells in which galaxies 
form and evolve. Over the past decades, observations have 
provided detailed information about the distribution of DM 
within those reg ions of s piral galaxies where th e baryons 



reside jAshmanl [19 92; Pers ic. Salucci & Stel 1 19961 (hereafter 
PSS); ISofue & Rubini.2001, ; , Salucci et a/.ll2007l) . Similar informa- 
tion on the distribution of DM is also ava i lable for low-surface- 
brigh tness (LSB) galaxies jde Blokl |2005| ; iKuzio de Narav et al\ 
I2OO6I) . In these disk systems, the ordered rotational motions and 
known geometry of the tracers has facilitated the mass modelling 
and provided clear evidence that the stellar components of spiral 
galaxies are embedded in extended DM haloes. In the most lumi- 
nous objects the stellar disk is almost self-gravitating with dark 



matter contributing significantly to the dynamics only at larger 
radii. In contrast, at the faint end of the galaxy luminosity function, 
baryons c ontribute a negligib le amount to the overall gravitational 
potential jde Blok et al 120081) . 



Extensive modelling of both individual and co-added spi- 
ral galaxy rotation curves (RCs) has generally concluded that 
almost maximal stellar disks embedded in cored dark matter 
haloes reproduce the data better than model s with cosmologically^ 



motivated, cusped dark matter haloes (PSS; [Persic. Salucci & Stel 
1996': 'Salucci & Burkert" 2000'; 'Gentile et a/."2004': 'de Blo k et al. 



2001 , 2002; Mai-chesin i 2002; Gentile 2005, 20d3, see 

also lcheminetalj|201ll) . Further, scaling relations between prop- 
erties of the spiral galaxies such as central surface density, stel- 
lar scale radius and stellar velocity dispersion have been identified 
and interpreted as signatures of the physical processes which drive 
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galaxy formation (e.g. | Kormend-vlll985l Il98l Il99d : lBurkertll99^ : 



iKomendv & Freemar 2004 ) 



iPersic & Salucci 



l l 19881) and lPersic & Saluccil ( ll99lh demon- 



strated that the Burkert halo density profile given by 



p{r): 



{r + ro) + rl) 



(1) 



with two free parameters, the core radius ro and the central halo 
density po, is consistent with the available kinematic data in spi- 
ral galaxies. When the mass distribution in t hese galaxies is mod- 
elled using the combined contributions of a iFreemanI ( Il970l) disk 
for the luminous matter and a Burkert profile for the dark matter 
halo, the structural parameters obtained (DM central densities, core 
radii, disk masses and length scales) exhibit a series of scaling laws. 
This led to the hypothesis of a "Universal Rotation Curve", an em- 
pirical function of radiu s and lumin osity that reproduces the RCs 
of sp i ral galaxies (PSS: jPersic & Salucci 1988 ; Persic & Saluccil 
Il99ll : lsalucci et a/.ll2007l and references therein). 

In contrast, our knowledge of the mass distribution in 
pressure- supported systems like elliptical galaxies is still lim- 
ited (see iNapolitano. Romanowskv & Tortora 1 1201(1 for a recent 
summary of the state of art). On-going observations of Local 
Group dwarf spheroidal galaxies (dSphs), which occupy the faint 
end of the luminosity function of pressure-supported systems, 
are currently yielding crucial information about the properties 
of the dark and luminous components in these objects and, in 
turn, on the underl yi ng physical p r opertie s of DM haloes (e.g . 
Tolstoy et al\ |2004|: iGilmore et al.\ I2OO7I: IStrigari e? a/ I l200p : 



Penarrubia. McConnachie & NavarrdhoOSi : IWalker et a/.H2009ah 



It is, however, an intrinsically difficult task both observationally, 
in terms of measuring velocities for sufficient numbers of tracers, 
as well as from a dynamical modelling point of view, due the lack 
of precise information on the dynamical state of the stellar popula- 
tions. 

The dSphs are typically at least two orders of magnitude less 
luminous than the faintest s pirals, and show ev idence of being DM 
dominated at all radii (e.g. iKlevna et al.ll2002l) . Their typical stel- 
lar masses lie in the range 3 x 10^ Mq to 2 x 10' Mq, although the 
lumino us masses of some re cently discovered objects are as low as 
10^ Mo jMartin et al.hooA while their stellar length scales are of 
order 0.3 kpc. In these systems, the DM halo typically outweighs 
the baryonic matter by a large factor (from a few tens, up to sev- 
eral hundred). An understanding of these objects is therefore es- 
sential for understanding the nature of dark matter itself and to 
build an observational picture of the outcome of galaxy forma- 
tion on small scales. Additionally, high-redshift dSph pre-cursors 
most likely contributed sig nifican tly to the build-up of the stel- 
lar halo of the Milky Way (,Helm i :2008 ^. Given that the observed 
dSphs are predominantly old, pressure-supported, spheroidal sys- 
tems, their evolutionary histories would be expected to differ sig- 
nificantly fro m those of spirals, e specially in the baryonic compo- 
nents (see e.g. lGrebel et al.l2003h . 

As in the case of spiral galaxies, a number of authors have 
found evi dence of univ ersality in the DM halo properties of dSph 
galaxies. iMated ( Il998h found that the variation of the mass-to- 
light ratios of dSphs with total luminosity was consistent with all 
dSphs containing similar masses of dark matter within the vol- 
ume occupied by their stellar distributions. This implies a larger 
proportion of dark matter in the less lu minous objects, a gen- 
eral characteristic o f spiral galaxy haloes l lPersic & Salucci 1 1 19881 : 
deBloket a"r]|2008l). More recent analyses jGilmore ef a/] |2007| : 



on extended velocity dispersion profiles rather than central velocity 
dispersions alone, have generally supported this conclusion. 

A number of important questions remain unanswered. First, is 
the distributio n of DM on galact ic (i.e. kpc) scales really universal? 
For example, lAden et al] ( l2009t) have noted the existence of con- 
siderable scatter in the estimated masses of the lowest luminosity 
systems, and several authors have presented evidence of systematic 
differences between the pro perties of the Milky Way dSph satellites 
and th ose surrounding M31 jCollins et alj201 lLl201ol : lKalirai et"al] 
l2010h . Secondly, why do the properties of the dark and luminous 
mass distributions appear to be related, even though baryons dom- 
inate, at most, only the inner regions of galaxies? 

The study of the internal kinematics of the Milky Way dSphs 
has been revolutionized by the availability of multi-object spec- 
trographs on 4m and 8m-class telescopes. Large data sets com- 
prising between several hundred and several thousand individ- 
ual stellar velocities per galaxy have now been acquired for all 



the luminous dSphs surrounding the Milky Wav JWilkinson et al. 
20041: iKlevna et fl/.ll2004l: iMunoz et fl/.ll2005l I2OO6I : IWalker et al. 



2009al : iKoch et a/.H2007il"bt iBattagha ef fl/.ll2008l) . The volume of 



the currently available data is sufficient to measure the dynamical 
masses interior to the stellar distributions of the dSphs. However, 
the mas s profiles remain les s well-determined. I t h as been demon- 
strate d l lWalker et al\ |2007|: iKoch et al\ l2007al lbl: iBattaglia et'al\ 
l2008h that these profiles are consistent with the cuspy dark mat- 
ter haloes produced in cosmo logical N-body simulations (e.g., 
iNavarro. Frenk & White Il997h , as well as with more general fam- 
ilies of haloes that r ange from centrally cored to steeply cusped 
jWalker et al .l2009bl) . The velocity dispersion profiles alone cannot 
distinguish between cored and cusped haloes du e to the degeneracy 
between mass and velocity anisotropy ( s ee e.g. iKoch et fl/.ll20p7al : 
iBattagli^ et a/.ll2008l : ISvans et al.|[2009l) . iGilmore et al.\ ( 120070 re- 
cently showed that the kinematic data and additional features in 
a small number of the well-studied dSphs are consistent with 
cored DM potentials, under the assumptions of spherical symme- 
try and velocity isotropy. In addition, several authors have pre- 
sented arguments which suggest that the internal kinematics of 



dSphs may be mo re consistent with cored h aloes ([Klevna et al 



Koch et a/.l2007al : IStrigari et al. l2008l : IWalker et Q/.l2009bl) , based 



I2OO2'; 'Goerdt et al."2006': 'Battag lia et ^l2008l : lAmorisco & Evans 
201 1 ; Walker & Pefiarrubia 201 ij). 

In the present paper, we investigate whether DM haloes of 
Burkert form are consistent with the observed kinematics of the 
luminous Milky Way dSphs. As noted above, Burkert haloes pro- 
vide good matches to the rotation curves of spiral galaxies and it 
is therefore interesting to ask whether they are also relevant mod- 
els for galaxies of other Hubble Types. We also wish to explore 
whether the parameters of the best-fit Burkert profiles for the dSphs 
lie on the extrapolation to the dSph regime of the scaling relations 
seen in spiral galaxies. While previous work has already shown that 
the internal kinem atics of the Milky W ay dSphs may be consistent 
with cored haloes jGilmore et a/.l2007l) . it does not necessarily fol- 
low that Burkert profiles in particular reproduce the observed kine- 
matic data. 

In what follows, we assume that the form of the dSph dark 
matter halo density is known and only the length scale and den- 
sity scale can vary. We allow the velocity anisotropy of the stellar 
distribution to vary in order to reproduce the observed dispersion 
profiles as closely as possible. Finally, we compare the resulting 
DM structural parameters with the low-luminosity extrapolation of 
the relations between the equivalent parameters found in spirals. 

Some of the comparisons between dSphs and spirals require 
us to define a stellar length scale for the dSphs which plays the 
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same role as the disk scale length Rq in spirals. One way to do this 
is to identify the location of the peak in the circular velocity curves 
in each system which would be predicted if the stellar components 
were assumed to contribute all the gravitating mass. For a Free- 
man disk, this peak occurs at 2 .2 R^. If the stel lar components of 
the dSphs are modelled using a |PlummejJl915l) sphere the corre- 
sponding radius is at 1.4 ri, (ri, is the projected half-light radius). 
Thus, where necessary, we associate the spiral disk length scale Ro 
with the radius 0.64 ri, in the dSphs. However, we note that most of 
our conclusions in this paper do not make use of this length scale. 

The outline of the paper is as follows. In Section|2l we summa- 
rize the observational data used in our study and describe in detail 
the analysis of the dSph data. Section [3] compares the properties 
of the dark haloes of spiral and dSph galaxies. Section |4] summa- 
rizes our findings and speculates on their implications for our un- 
derstanding of dark matter and galaxy formation. 



2 DATA 

2.1 dSph galaxies 

Figure [T] displays e mpirical velo c ity di s persion profil e s orig - 
inally publis h ed by IWalker et all ll2007h: iMateo et al. I feOOsh : 
I Walker ef al\ l l2009ah and IWalker et al] ( l2009bl) for the Milky 
Way's "classical" dSph satellites Carina, Draco, Fornax, Leo I, Leo 
II, Sculptor, Sextans and Ursa Minor. As discussed above, in our 
present analysis, we assume Burkert profiles for the dark matter 
haloes of the dSphs in order to provide a basis for comparison with 
spiral galaxies. Specifically, we use the previously published veloc- 
ity dispersion profiles shown in Figure [T] to constrain values of the 
Burkert parameters po and ro for each dSph in our sample. 

We assume that the luminous component of each dSph con- 
sists of a single, pressure-supported stellar population that is in 
dynamical equilibrium and therefore traces the underlying gravi- 
tational potential which we assume to be dominated by the dark 
matter halo. The masses of these stellar spheroids can be estimated 
from their luminosities: they are 1-2 orders of magnitude smaller 
than the dynamical masses. We assume that the stellar mass-to-light 
ratios (M/L)v are unity. While the actual (M/L)v ratios may vary by 
about a factor of two depending on the details of the stellar popula- 
tions (see e.g. lMateo et al .119981) , the uncertainties in the modelling 
results are dominated by the unknown velocity anisotropy. 

The Jeans equation relates the density and velocity disper- 
sion of the stellar component to the mass profile of the dark matter 
halo. Assuming spherical sym metry, the Jeans equation f or a non- 
rotating system can be written teinnev & TremaiMl2008l) 



V ar r 



GM{r) 



(2) 



where v(r), v;(r), and j8(r) = 1 - v^/i'^ represent the 3-dimensional 
density, radial velocity dispersion, and orbital anisotropy, respec- 
tively, of the stellar component, and M(j) is the mass profile of 
the dark matter halo. In this analysis the orbital anisotropy fUj) is 
not constrained, as all information about the velocity distribution 
is restricted to the component along the line of sight. We make the 
simplifying assumption that ji = constant, which prov ides the fol- 
lowing solution to EguationEll lMamon & Lokas II200I : 



v{s)M{s)ds. 



(3) 



In order to compare to observable s, we consider the project ion of 
Equation[3]along the line of sight l lBinnev & Tremaingl2008l) : 



j_ r 

m) Jr 



R^ 1 v\'fr 



) Vr2 - R^- 



(4) 



where I{R) is the projected stellar density profile and (T,,(^) is the 
projected velocity dispersion profile. The two parameters of inter- 
est are of course the central density and core radius, which enter 
Equation |4] upon substituting for vv^. (Equation [3]l with the mass 
profile derived from the Burkert density profile: 

M{r) = 4n [ s-p(s)ds 
Jo 

= npori (ln[(I + r/rofd + r^/r^)] - 2tan-'[r/ro]) . (5) 

To describe the stellar density profile we adopt a Plummer 
profile, I(R) = L(7:r^)-^[l + R-jrlY^, where L is the total lu- 
minosity and ri, is the projected half-light radius (i.e., the radius 
of the circle that encloses half of the total luminosity in projec- 
tion). Under the assumption of spherical symmetry, the correspond- 
ing 3-D st ellar density i s then v(r) = 3L(4;rrj^)"'[l -1- r-jr^Y^^^-. 
Following IWalker et al] ( l2009bl) . for the eight dSphs considered 
here we adopt the V-b a nd lum inosities and half-light radii from 
llrwin & Hatzidimitrioul ( [l"995h . All values are tabulated in Table 
I of lWalkere?fl/.l ( l20I(]|) . 

Treating the stellar density as a known function, we fit 
halo models to the empirical velocity dispersion profiles us- 
ing the set of three free parameters: 9 = [61,62,63] = 
|logio[''o/pc],log,o[po/(Mopc-3)],-logio(l -yS)). We adopt uni- 
form priors over the ranges -10 < log(()[po/(Mopc"^^)] < 5, 
-2 < log[o[ro/pc] < 5, and -1 < -log,o(l -y8) < 1. For a given 
point in parameter space. Equation |4] specifies the projected veloc- 
ity dispersion profile crp{R). We compare model profiles to the em- 
pirical velocity dispersion profiles, crvf,{R) (FigureQ}, by evaluating 
the likelihood 



L(0) = Y] 



: exp 



]^(<tI(R.) - o-liRdf 

'2 



Var[a-2^(i?,)] 



(6) 



where Var[(Ty^ (/?,)] is the square of the error associated with the 
square of the empirical dispersion and N is the number of bins 
in the dispersion profile. We obtain (marginalised) 1-D posterior 
probability distribution functions for each free parameter using 
a Markov-Chain Monte Carlo (MCMC) algorit hm. Specifically, 
we use the same Met ropolis-Hastings algor ithm jMetropolis et ali 
ll953l : lHastin"gs]|l970l) described in detail bv lWalker et al] l l2009bh . 
In order to account for the error associated with observational un- 
certainty in the half-light radius, for each point sampled in our 
MCMC chains we scatter the adopted value of by a random de- 
viate drawn from a Gaussian distribution with standard dev iation 
equal to the published error jlrwin & Hatzidimitrioul [l99^ . The 
probability distribution functions for the free parameters are thus 
elfectively marginalised over the range of half-light radii consistent 
with observations. 

The velocity dispersion profiles corresponding to the highest- 
likelihood (Equation |6j point from each of our MCMC chains are 
over-plotted on the empirical profiles in Figure [T] These 'best fits' 
demonstrate that Burkert profiles can provide an excellent descrip- 
tion of dSph velocity dispersion profiles. For each free parameter, 
we take the 1-D posterior probability distribution obtained from our 
MCMC chains as the observational constraint, given our modelling 
assumptions. For each free parameter (and combinations thereof), 
Table[T]identifies the median value and confidence intervals that en- 
close the central 68% (and 95%) of accepted points in our chains. 
We find that the dSph haloes have central densities ranging from 



4 P. Salucci et al. 



15 

10 

5 


15 

10 

5 



' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' 
-Burkert Carina 




Draco - 



_L 



250 



500 



750 



1000 



500 



1000 



1500 



2000 



_ 11 11 1 


111,111 


, 1 , , , 1 _ 
Fornax - 












— : 


- , , , , 1 




, 1 , , , , - 



Leo 




1000 



500 



250 500 750 1000 1250 
R [pc] 



250 500 750 1000 
R [pc] 



Figure 1. Velocity dispersion profiles for the Milky Way's eight "classical" dSph satellites. Over-plotted are the best-fitting profiles obtained under the 
assumptions of Burkert dark matter haloes, Plummer light profiles and radially constant velocity anisotropy. The parameters of each fit, together with associated 
confidence limits, are listed in Table[T] 



cm"^ to 3X lU""^g cm"'' and core radii ranging from 0.05 
kpc to 0.65 kpc. The data in the table also exhibit the well-known 
mass-anisotropy degeneracy: because the dispersion profiles of the 
dSphs are essentially flat to large radii and have a relatively small 
range of amplitudes, larger values of po''o are required for galax- 
ies with more radially biased velocity distributions. Our analysis is 
particularly susceptible to this degeneracy due to our restriction of 
the modelling to Burkert halo profiles. 



We have repeated our analysis with the additional assumption 
of velocity isotropy (i.e. [} = 0). With the exception of Sextans, 
the halo parameters obtained for our sample are consistent with 
those in Table [T] within the quoted uncertainties. When we restrict 
ourselves to isotropic models, the best-fit rg for Sextans falls to the 
significantly smaller value of 47pc. This is consistent with the fact 
that Sextans is unique in requiring tangential anisotropy to obtain a 
good fit to the dispersion profile - the best-fit isotropic model does 
not match the profile of Sextans interior to 200pc. 



2.2 Spiral Galaxies 

As discussed above, Burkert halo models provide excellent fits 
to individual spiral galaxy rotation curves as well as to sam- 
ples of co-added rotation curves. Moreover, when the mass mod- 
elling is performed using Burkert haloes, a tight relation between 
po and To (and also between ot her parameters like the disk and 
virial masses) eme rges (PSS; iDonato. Gentile & Saluccil |2004 
ISalucci et a7ll2007h . As can be seen in Figure |2] we find similar 
Po vs ro relationships independently of whether the mass profiles 
are obtained from rotation curves or from gravitational lensing data 
and irrespective of whether the analysis is performed on individual 
or co-added rotation curves. 

To emphasise the very different ranges of baryonic mass and 
extent probed by the dSphs and the spiral galaxies in our sample, 
in Figure[3] we compare the relationship between the characteristic 
baryonic length scale (Rj^; see above for definitions) and the stellar 
mass of dSphs (estimated from the V band luminosity, assuming a 
stellar mass-to-light ratio of unity [in solar units]), and of spirals. 
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Table 1. Tabl e of dSph structural parameters and results of our mass modelling. The columns are: (1-3) dSph name, observed V-band luminosity and half-light 
radius (from llrwin & Hatzidimitrioj[l99^ ; (4) velocity anisotropy; (5-6) density scale and length scale for the DM halo, assuming a Burkert DM density 
profile; (7) halo "surface density" scale; (8) total mass interior to Rg3/2, where Rg^ is the radius enclosing 0.83 of the total luminosity; (9) ratio of stellar mass 
to DM mass interior to Rg3/2. In columns (4-9) the median value of the parameter is given together with error bars which enclose the central 68% (95%) of the 
marginalised ID probability distribution function for the pai'ameter. Since we assume that M/Lv=l for the stellar components of the dSphs, the luminosities 
directly yield the stellar masses. 
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Figure 2. Pai'ameters of Burkert DM haloes obtained from dynamical mod- 
elling of (i) spirals, based on the URC hypothesis applied to co-added ro- 
tation curves (solid line; data from PSS) or weak lensing shear (squares; 
iHoekstra et"aDl2005l) ; (ii) NGC 3741 (triangle) the darkest spiral in the lo- 
cal Universe based on its kinematics; (iii) the "classical" Milky Way dSph 
satellites (filled circle s), based on their internal stellar kinematics (this pa- 
per). The lSpano etal\ 1200 8i) relation is shown as a dashed line. 

3 DARK MATTER SCALING RELATIONS 

In the previous section, we showed that Burkert halo profiles pro- 
vide good fits to the dSph kinematic data (Figure[T]l. In this section 
we compare the parameters of the Burkert profiles obtained from 
our dynamical modelling of dSphs with those obtained for spiral 
galaxies. In Figure (2] we plot po versus ro for the eight dSphs in 
our sample. Remarkably, they lie on the extrapolation to higher 
central densities of the relation found for spirals. All these data 
can be reproduced by the relation log po - log ro +const with 
0.9 < < 1.1. 

An even more interesting comparison can be done involv- 
ing the mean dark matter surface density within the dark halo 
core radius (the radius within which the volume density profile 



of da rk matter remains ap proximately flat). It was recently discov- 
ered jPonato et al\\200<^ that this quantity fiQ (jj.q = paro) is con- 
stant across a wide range of galaxies of different Hubble Type and 
luminosity and that this relation holds also for dSphs, as we confirm 
in Figure|4]where the data for the dSphs ar e the results of o ur Burk- 
ert halo modelling in the present paper (in lDonato et al. I the dSph 
halo parameters were obtained via a different approach). We there- 
fore confirm that this relation extends across a luminosity range of 
14 magnitudes and spans the whole Hubble sequence . The potential 
implications of the constancy of /jq are discussed in lDonato et al. I 
(2009). 

In their modelling of spiral galaxies using the URC hypothe- 
sis, PSS found that the parameters of their Burkert DM haloes were 
correlated with those of the luminous matter. In Figure [5] we show 
the Po vs Ro relationship for our dSph sample compared with the 
corresponding relation for spirals from PSS. As in the case of the 
po-''o relation in Figure|2l the dSph data are consistent with the ex- 
trapolation of the relation seen in spirals. The significance of this 
relation derives from the fact that it links the DM and baryonic mat- 
ter properties of galaxies on a wide range of length scales: qualita- 
tively, the "central" densities of DM haloes increase as the extents 
of their associated stellar components decrease Although the obser- 
vational evidence for this relation is relatively strong, we stress that 
its physical interpr etation is presently unknown (see I Gentile et al.l 
l2009l ; lAngusll2008l. for some related discussion of this point). 

It is interesting to consider explicitly the role played by veloc- 
ity anisotropy in this result. Unsurprisingly, under the assumption 
of Burkert haloes for dSphs the inclusion of velocity anisotropy as 
a free parameter improves the quality of the dispersion profile fits 
relative to those obtained for isotropic models. However, we also 
find that the scatter in the po-ro relation is smaller for anisotropic 
models - thus the better we reproduce the observed dispersion pro- 
files using Burkert haloes, the tighter is the correlation between the 
halo parameters. 

Finally, we emphasise that the present results, do not require 
the presence of cored haloes in dSphs, nor do they constrain the 
density and scale lengths of their haloes in a model-independent 
way. On the other hand, the fact that the dSph kinematics can be 
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Figure 3. Comparison of the distribution of characteristic baryonic scale i?D 
versus stellar mass Ms for dSphs (points; this paper) with the corresponding 
relation in Spirals (from PSS). See Section[T]for the definition of Sd used 
for the dSph sample. 




ferent galaxies and Hubble types. The data are: (1) the lspano et al\ lloOSh 
sample o f spiral galaxy data ( open red circles); (2) the URC relation (solid 
blue line:IShank:ar et all2006l) ; (3) the dwarf in'e gulars N 3741 (Mb = 13.1; 
iGentile et al II2007I) and DDO 47 (Mb = 14.6; iGentile et al II2005I) (ful l 
green circles), spirals and ellipticals (black squares; iHoekstra et al.ll2005i) 
investigated by weak lensing; (4) Milky Way dSphs (pink tri angles - this 
paper ); (5) nearby spirals in THINGS (small blue tri angles; IWalter et alj 
20081); (6) early-type sp kals (full red triangles; lNoorderm eei[_[2006l; 
Noordermeer et alTl2007l) . The long-dashed line shows the iDonatoeTa/n 
<200g|) result. 



reproduced using Burkert DM halo profiles whose structural pa- 
rameters lie on the same scaling relations as those of spirals pro- 
vides some support for the claim that the mass distributions in dSph 
galaxies can be understood within the same framework as those of 
spirals. 



4 CONCLUSIONS 

Dwarf spheroidal galaxies are the lowest luminosity stellar systems 
which show evidence of dynamically significant DM. Their physi- 
cal properties (luminosity, stellar scale length, baryon fraction) are 
typically two orders of magnitude different from those observed for 
spiral and elliptical galaxies. Given these extreme structural prop- 
erties, an understanding of the formation of dSphs is crucial for the 
development of a complete picture of galaxy formation. 



o " 
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Figure 5. Halo central density po versus stellar length scale for spirals 
(solid curve) and dSphs (points). 



The main result of this paper is the finding that these galaxies, 
despite being very distinct in their physical properties from spirals 
and ellipticals and having a large individual scatter in their baryonic 
properties, exhibit kinematic properties that can be modelled using 
DM haloes with the same mass profiles as those which reproduce 
the rotation curves of spiral galaxies. Under the assumption that 
the haloes of dSphs have Burkert profiles, we find that the derived 
central densities and the core radii are consistent with the extrap- 
olation of the relationship between these quantities seen in spiral 
galaxies. Conversely a Burkert profile with structural parameters 
predicted by the extrapolation of the relation between halo central 
density and DM core radius previously found from Burkert fits to 
the kinematics of elliptical and spiral galaxies can account for the 
observed internal kinematics in dSphs. 

This result is intriguing, and could point to a common phys- 
ical process responsible for the formation of cores in galactic 
haloes of all sizes, or to a strong coupling between the DM 
and luminous matter in dSphs. It is worth noting that a poten- 
tial connection between spiral galaxies and dSphs does not ap- 
pear as natural as one between dSphs and other hot, s pheroidal 
systems jPabringhausen etal. 1 l2008l ;l Forbes et al. I I2OO8I) . For ex- 
ample, while the sizes of spiral galaxies are presumably fixed by 
the angular momentum of the gas from which they form, most 
of the present -day dSphs show no signs of rotation (although 
iBattaglia et al\ 1 120081) have r ecently found evidence of rotation in 
the Sculptor dSph). However. iMaver et al. I ( |20()]1) have proposed a 
formation scenario for dSphs in which they are initially low-mass 
disk galaxies that are subsequently transformed into spheroids by 
tidal interaction with the Milky Way. More recently, such models 
have been sho wn to provide reasonable m odels for t he properties 
of the Fornax jKlimentowski et a/ni2007h and Leol fcokas et al. I 
l2008h dSphs. If the haloes of dSphs do indeed follow the scaling 
laws defined by more massive disk galaxies, this could lend indirect 
support to evolutionary histories of this kind. Suggestive evidence 
of such transformation scenarios is also provided by the discovery 
of residual disks with spiral structure in lumin ous dwarf elliptical 
galaxies in the Virgo cluster jUisker et al.l2006l) . 

Further dynamical analysis is needed to derive the actual DM 
distribution in dSph and possibly to estimate their halo core radii. 
Nevertheless, it interesting to speculate on the possible implications 
of these scaling laws for our understanding of DM. Warm dark mat- 
ter has been invoked as a potential solution to the over-prediction of 
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substructure by AC DM simulations, and to tiie cusp-core issue (e.g. 
iMoore et al. 1 19991) . However, the existence of scaling relations be- 
tween the central density and core radius over three orders of mag- 
nitude in both quantities would argue against this explanation, un- 
less the warm DM spectrum is extremely fine-tuned. Further, such 
DM relations cannot arise due to either self-annihilation or decay 
of DM which would predict a narrow range in po and no clear cor- 
relat ion of the latter with t he cor e radius. 

iDalcanton & Hogaiil ( 1200 ih argued that the phase-space den- 
sities of DM haloes suggested that warm DM (either collisional or 
coUisionless) could not be the cause of cores in galaxy haloes on 
all scales. These authors suggested a dynamical origin for the cores 
of larger galaxies. Subsequently, a number of studies have demon- 
strated that macroscopic core formation in galaxy haloes is possible 
through the infall of compact baryonic (or baryon dominated) sub- 
clumps (El-Za nt et alJlOOlUjardel & Sellwoodll2009l ; lGoerdt et al.l 
lioTo iCole et al.l201 l|). Further work is required to explore whether 
such processes , in conjunction wi th subsequent star formation and 
feedback (e.g. IPasetto et al.ll2010l) . can result in universal scaling 
relations spanning three orders of magnitude in density and length 
scales. 

Clearly, a direct kinematic determination of the dark mat- 
ter profiles of dSphs is essential to confirm the robustness of the 
scaling relations between the halo parameters. A number of re- 
cent papers have made progress towards this goal in the subset 
of dSphs which exhibit kinematically distinct sub-population s in 
their stellar components tBattaglia e t al. 2008; Amorisco & Evans! 
I2OI ll : I Walker & Penarrubiall201ll) . Interestingly, all three analyses 
favour cored haloes over cusped ones. However, further observa- 
tional and modelling work is required to constrain the halo profiles 
more tightly. Theoretically, the development of a physical picture 
of the processes which shape the halo profiles of dSphs and which 
could lead to the existence of apparently similar scaling relations 
between halo properties over a wide range of galaxy luminosities, 
is an important research goal for the coming years. 
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